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Context. Tori of Active Galactic Nuclei (AGN) are made up of a mixture of hot and cold gas, as well as dust. In order to protect the 
dust grains from destruction by the surrounding hot gas as well as by the energetic (UV/optical) radiation from the accretion disk, the 
Q\ ' dust is often assumed to be distributed in clouds. 

Aims. A new three-dimensional model of AGN dust tori is extensively investigated. The torus is modelled as a wedge-shaped disk 
within which dusty clouds are randomly distributed throughout the volume, by taking the dust density distribution of the corresponding 
continuous model into account. We especially concentrate on the differences between clumpy and continuous models in terms of the 
temperature distributions, the surface brightness distributions and interferometric visibilities, as well as spectral energy distributions. 
Methods. Radiative transfer calculations with the help of the three-dimensional Monte Carlo radiative transfer code MC3D are used 
q , in order to simulate spectral energy distributions as well as surface brightness distributions at various wavelengths. In a second step, 

j_( . interferometric visibilities for various inclination as well as position angles and baselines are calculated, which can be used to directly 

' compare our models to interferometric observations with the MIDI instrument. 

' Results. We find that the radial temperature distributions of clumpy models possess significantly enhanced scatter compared to the 

continuous cases. Even at large distances, clouds can be heated directly by the central accretion disk. The existence of the silicate 
10 |rm-feature in absorption or in emission depends sensitively on the distribution, the size and optical depth of clouds in the innermost 
part of the dust distribution. With this explanation, failure and success of previous modelling efforts of clumpy tori can be understood. 
The main reason for this outcome are shadowing effects of clouds within the central region. We underline this result with the help 
of several parameter variations. After adapting the parameters of our clumpy standard model to the circumstances of the Seyfert 2 
Circinus galaxy, it can qualitatively explain recent mid-infrared interferometric observations performed with MIDI, as well as high 
• resolution spectral data. 
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1 . Introduction and motivation is blocked and optical emission lines can only be detected from 

gas beyond the torus funnel. Being further away from the centre, 

According to today s knowledge, Active Galactic Nuclei (AGN) ft moyes slower and hence duces narrow emission lines onl 

are powered b y accretion ont o a^upermassive black hole (10 6 - This is the so-ca l led Unified Scheme of Active Galactic Nuclei 

H . 10 10 M , e.g. I ShankaretalJ |2004j) residing in their centres. dAntonuccill993blUrrv & Padovanil 19951) . First evidence for this 

Thereby, gravitational energy is converted into heat by viscous came from S p ec trop olarimetric observations of type 2 

processes within the surrounding accretion disk, which extends sources ,|Miller & Antonucd fl98lh . clearly displaying type 1 

from the marginally stable orbit up to several thousands of ^natures in the polarised light, which is scattered by electrons 

Schwarzschild radii. The emitted UV/optical light illuminates and tenuouf , dust wkhin ^ funnd aboye ±e tQms The opening 

the attached, toroidally shaped dust reservoir. The concept of le of the toms can be estimated w ith the help of stati stics of 

this obscuring torus was introduced in order to unify main y two ±e different of § fert laxief , Maiolino & Riekel (fl995l) 

classes of observed spectral energy distributions (SEDs): one find a ratio between s 2 to s x laxies of 4:1 in their sam . 

shows a peak in the UV-region with overlayed broad and narrow ^ which resuhs fa an in k of ^ R U CQnes of ^ 

optical emission lines, the other class shows only narrow optical ^ concordance with observations of ionisation cones of 

emission lines. This can be interpreted as an inclination angle individual alaxies Direct support for the idea of geometrically 

dependence. For viewing angles within the dust-free cone of the Mck torf CQmes f romrecent interf erometric obser vations in the 

toms (type 1 sources) direct signatures of the accretion disk (a mid . infrared (e |j affe et alf2 004: Tri stram et al.ll2007h . 
peak in the UV-range) and the region close to the centre within 

the funnel of the torus show up. This is where gas moves fast and, These tori are made up of at least three components: (i) hot 

therefore, produces broad emission lines (the region is hence ionised gas, (ii) warm molecular gas and (iii) dust. Krolik & 

called the Broad Line Region (BLR) of the nucleus). For edge- Begelman (1988) proposed that the dusty part has to be organ- 

on lines of sight (type 2 sources), the direct view onto the centre ised in a clumpy structure in order to prevent the grains from 

being destroyed by the hot surrounding gas (with temperatures 

Send offprint requests to: M. Schartmann of the order of 10 6 K) in which the clouds are supposed to be 

* e-mail: schartmann@mpe.mpg.de embedded. Another hint for the clumpy nature of the obscuring 
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material - in this case mainly for the distribution of neutral gas - 
com es from X-ray measu rements of the absorbing column den- 
sity. iRisaliti et al.1 (|2002) claim that the observed variability of 
these measurements on timescales from months to several years 
can be explained by a clumpy structure of the torus. Combining 
X-ray absorbing column densities with spectral information fur- 
ther strengthens the claim for a clumpy distribution of the dust 
(IShi et alJl2006h . 

Earlier work on torus simulation s concentrated mostly 
on smooth dust distrib ut ions (e.g. Pier & Krolikl [l992a; 
Granato & Danesd 1 1994 Ivan Bemmel & Dullemondl 120031: 



Schar tmann et al. 1 120051) . This was mainly caused by the lack 



of appropriate (3D) radiative transfer codes and computational 
power. Nevertheless, such models are good approximations 
for the case that the clumps that build up the torus are small 
compared to the total torus size, as is also shown in a parameter 
study described in Sect. 14.21 These continuous models are 
able to de scribe the gross observa ble features of these objects 
(see e. g. ISchartmann et al.l 120051) . However, problems arose 
from too strong emission features of silicate dust compared 
to the observations, when looking directly onto the inner rims 
of the model structures (face-on views). They had never been 
observed before that time, although almost all models showed 
them for the face-on view. Therefore, much theoretical effort 
was undertaken in order to find models showing no silicate 
feature at all in the face-on case, while retaining the s ilicate 
absorption feature in the edge-on case. iManske et al.l dl998l) 
for example succeeded in avoiding silicate emission features 
with a flared dust disk of high optical depth in combination 
with an anisotropic radiation characteristic of the central 
illuminating source. A very promising idea was to solve the 
problem naturally by splitting the d ust distribution i nto sin gle 
clouds. This was first attempted by iNenkova et al.l (12002b . A 
one-dimensional code for the simulation of radiative transfer 
through single clumps was used and, in a second step, the 
torus and its emitted SED was assembled by adding many 
clouds of different aspect angles with the help of a statistical 
method. With this approach, they could show that a clumpy 
dust distribution of this kind can significantly smear out the 
prominent silicate emission feature of the SEDs of type 1 
objects at 10 \im for a large range of parameter values. No more 
fine-tuning was needed, as in the previously proposed solutions 
with the help of special continuous models. Subsequently, real 
tw o-dimensional radiative transfer ca lculations were undertaken 
bv lDullemond & van Bemme]| ( l2005l) . Clouds were modelled as 
concentric rings. A direct comparison between these kinds of 
clumpy models and the corresponding continuous models did 
not show evidence for a systematic suppression of the silicate 
feature in emission in the clumpy models. 

Meanwhile, silicate features in emission were found with 
the help of the Infrar ed Spectrograph (IRS) onboard the Spitzer 
space telescope (e . g.lSie benmorgen et al. 2 0051 lHao et al.l2005t 
ISturm et al.l 120051: IWeedman et all 120051) . For these kinds of 
studies, Spitzer is superior to other available facilities, due to its 
high sensitivity and the coverage of a wavelength range includ- 
ing both silicate features (at 9.7 um and 18.5 um) and the sur- 
rounding continuum emission. Silicate emission features were 
found in different levels of AGN activity, ranging from very lu- 
minous quasars down to weak LINERS. These findings are in 
good agreement with a geometrical unification by an optically 
thick dusty torus, as silicate emission features can be produced 
even in the simplest models. But one has to be cautious, as due 
to the large beam of the Spitzer space telescope and the low tem- 
peratures measured, it is unclear whether these silicate features 



result from dust emission in the innermost parts of the torus or 
from optically thin regions surrounding them. 

Very detailed simu l ations of clumpy tori were undertaken re- 
cently by Honig et al.l (120061) . They apply a similar method as 
INenkova et al.l (120021) . but use a 2D radiative transfer code for 
the simulation of SEDs of individual spherical clumps at var- 
ious positions in the torus and with various illumination pat- 
terns: directly illuminated and/or illuminated by reemitted light 
of surrounding clouds. In a second step, these clouds are dis- 
trib uted according to physica l models by IVollmer et al.l (120041) 
and iBeckert & Duschll ( 120041) . A comparison of the resulting 
SEDs and images with spectroscopic and interferometric obser- 
vations shows good agreement. This model is characterised by a 
large number of small clouds with a very large optical depth, es- 
pecially close to the centre. We compare our models with these 
models in Sect. 15. 21 

Despite the detection of geometrically thi c k dust tori in 
nearb y Seyfert galaxies (e.g. iJaffe et all 120041: iTristram et al.l 
120071) . many questions remain: How are these tori formed? How 
are they stabilised against gravity? Do steady torus solutions ex- 
ist? Several attempts to answer th e se que stio ns have been made. 
For ex ample iKrolik & Begelmanl dl988l) and IBeckert & Duschll 
(|2004|) support the scale-height of their tori with the help of 
discrete clumps, moving at supersonic velocities, maintained 
by mainly elastic collisions with the help of strong magnetic 
fields. Other gr oups replace the torus b y a magnetically-driven 
wind solution (iKonigl & Kartid 1 19941) . The most recent sug- 
gestion comes from jKrolikl (12007b . building up on an idea of 
iPier & Krolikl (fl992bh~ where the scale-height of tori can be 
maintained with the help of infrared radiation pressure, as shown 
with an idealised analytical model. A more det ailed review o f 
possible solutions and their drawbacks is given in Krolik 120071) . 
Another possible scenario, where the effects of stellar feed- 
back from a nuclear clu ster play a major role, is discussed in 
ISchartmann et al l d2008l) . 

In this paper, we address the implications of dumpiness on 
the temperature structure, the infrared spectral energy distribu- 
tions, surface brightness distributions as well as interferometric 
visibilities by implementing fully three-dimensional radiative 
transfer calculations through a clumpy dust distribution and 
discuss the possible mechanisms causing this behaviour. In 
Sect. [2] a description of our model is given, before we present 
the basic results for our standard model (Sect. [3} and for several 
parameter studies (Sect.|Ui and discuss the findings (Sect.|5]l, 
as well as differences and similarities to other models. In 
Sect.[6]we interpret our results in terms of MIDI interferometric 
observations and compare them to data for the Circinus galaxy. 
Finally we draw our conclusions in Sect. [7] 



2. The model 

2. 1 . Assembly of our clumpy standard model 

We apply a very simple, wedge-like torus geometry with a half 
opening angle of 45° in order to gain resoluti on. In our previous 
two-d imensional continuous 7TM-models ( ISchartmann et al.l 
I2005I) . the simulation of the whole 0-range was necessary, due 
to the radial as well as 0-dependence of the dust distribution. It 
resulted from an equilibrium between turbulent pressure forces 
and forces due to an effective potential. The latter is mainly 
made up of gravitational forces due to the central black hole and 
the central stellar distribution, as well as rotation. The cloudy 
dust distribution is set up on a spherical three-dimensional grid 
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Table 1. Main model parameters for our continuous and clumpy 
standard model. 



both models 



inner radius of the torus 

outer radius of the torus 

half opening angle of the torus 

total optical depth in equatorial plane 

exponent of continuous density distribution 
number of grid cells in r direction 
number of grid cells in 6 direction 
number of grid cells in <p direction 



^out 
Vat 



L'open 
eqi 



0.4 pc 
50 pc 
45° 
2.0 

-0.5 
97 
31 
120 



additional in clumpy model 



number of clumps 
exponent of clump size distribution 
constant of clump size distribution 
optical depth of each clump 
average number of cells per clump 



P 

a 

clump 
T 9.7pm 



400 
1.0 

0.2 pc 
0.38 
272 



r = (r, 9, (f>), which is linear in 9 and <p an d logarithmic in r. 
To obtain the clumpy density structure, the following proce- 
dure is applied: A random number generator (RAN2 taken from 
iPress et alii 19921) determines the radial coordinate of the clump 
centre, which is equally distributed between the inner and outer 
radius. The 9 and (p coordinates are chosen such that the result- 
ing points are equally distributed on spherical shells. In a second 
step, the spatial distribution found so far is coupled to the dust 
density distribution of the continuous model: 



(1) 



Pcont(r, 9, cp) = po\ 

,lpc 



The radii of individual clumps a c i ump vary with distance from 
the centre according to a distribution 



^clump 



(^"clump \ 



(2) 



All cells within this clump radius are homogeneously filled with 
dust. All clumps possess the same optical depth t^' 1 ™^, measured 
along a radial ray through the clump centre. We further require 
that clumps always have to be completely contained within the 
model space, but are allowed to intersect. Such a combination of 
intersecting clumps will be called a cloud from now on. Thus, a 
cloud may contain overdensities, where the intersection happens. 
A clump size distribution as described above seems to be rea- 
sonable, as shear forces due to the differential rotation increase 
towards the centre. Therefore, clouds are more easily disrupted 
in the inner part of the torus. Furthermore, clouds become com- 
pressed when moving towards the centre due to the increasing 
ambient pressure in a deeper potential well. 

All other routines and algorith ms used in this paper are i den- 
tical to the modelling described in Schartmann et al. (2005) and 
will only be mentioned briefly in Sect. |2.2| 

The main model parameters of the continuous and clumpy 
distributions are summarised in TableQ] where the numerical 
values refer to our clumpy and continuous standard model. The 
dust density distribution for the clumpy case is shown in Fig.[T] 
The torus possesses a volume filling factor of 30% and the 
dust mass was chosen such that the optical depth of the torus 
within the equatorial plane (averaged over all angles <p) reaches 



Fig. 1. 3D rendering of the clump distribution of our standard 
clumpy torus model. The chosen inclination angle corresponds 
to a Seyfert 1 type (face-on) view onto the torus. 



a value of two at 9.7 um. With this value, the resulting absorp- 
tion column densities are in concordance with observations of 
Seyfert type 2 galaxie s obtained with the IRS spectrometer on- 
board Spitzer (see e. g. lShi et al.l2006l) . and the modelled silicate 
absorption feature depth compares well with observations. If not 
stated otherwise, the optical depth always refers to a wavelength 
of 9.7 u rn throughout this pape r. 

As in lNenkova et al.l d2002l) . all clumps possess the same op- 
tical depth in our standard model. A total optical depth within the 
equatorial plane of (t^^ = 2.0 results in an optical depth of 

T 97™m = 0-38 along a radial ray through the centre of the clump. 
The corresponding continuous model has the same geometrical 
structure, continuously filled with dust according to the density 
distribution given in equationQ] 



2.2. Summary of methods 

A very brief overview of the dust composition, the heating 
source and the numerical method of radiative transfer will be 

given in this section. 

Although several hints (iMaiolino et al.l 12001 alibi iJaffe et al.l 
|2004|) point to the possibility that dust in the nuclear regions of 
AGN is dominated by large grains, we will limit our presen t in- 
vestigation to the classic MRN-model dMathis et all 1 19771) for 
three reasons: first, and most important, we ai m for comparabil- 
ity w ith our earlier paper on continuous tori dSchartmann et al.l 
120051) . Second, we have tested that our essenti al results about the 
chang e in grain size distribution (Sect. 3.9 in ISchartmann et aTl 
120051) remain unchanged when distributing the dust in a clumpy 
structure. Third, our approach, which explicitly takes into ac- 
count the size-dependent sublimation radius, is generically more 
robust against changes in the grain distribution than calculations 
that ignore this effect. For our current simulations, we represent 
the MRN-model by three different grain species with 5 different 
grain sizes each. Taking different sublimation radii of the various 
grains into account then partially accounts for the destruction of 
small grains in the harsh environment of the quasar, as they pos- 
sess larger sublimation radii. 
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Fig. 2. a) SEDs for a photon number study. The solid curves re- 
fer to our standard model and the dotted graphs (identical with 
the solid curves) result after doubling the number of photon 
packages, b) SEDs for a resolution study: high resolution (solid 
curves - our standard model) and a factor of 3 reduced number 
of grid cells (dotted curves). Shown are the cases for inclination 
angles 0° and 90°. 



The dust distribution is heated by a point-like, central accre- 
ti on disk with the SEP of a mean quasar spectrum (see Fig. 3b 
in lSchartmann e t al. 20050. The radiation characteristic is chosen 
to follow a | cos(0)| law for all wavelengths. For the simulations 
shown in this paper, the accretion disk SED is normalised to a 
bolometric luminosity of 1.2 x 10 11 Lq, except for the compari- 
son with the Circinus galaxy. 

In order to obtain the temperature, the SEDs and the sur- 
face brightness distributions of the dusty torus, we use the three- 
dimensiona l radiative transfer code MC3EQ] ( I Wolf et alj|1999t 
IWolfll200"3l) . We apply the Monte Carlo procedure mainly for the 
calculation of temperature distributions and the scattering part 
whenever necessary, whereas SEDs and surface brightness maps 
for dust reemission are obtained with the included raytracer. The 
main advantage compared to other codes is MC3D's capability 
to cope with real three-dimensional dust density distributions, 
needed for a realistic modelling of the dust reemission from 
a clumpy torus. For this paper, we implemented the automatic 
determination of the sublimation surfaces of the various grain 
species in three dimensions. As we expect the sublimation to 
happen along irregularly shaped surfaces in a three dimensional, 
discontinuous model, a raytracing technique is used to solve the 
(ID) radiative transfer equation approximatively in all directions 
of the model space. 

For further information on the radiative transfer proce- 
dure used and the other pr econditions (main l y primary source 
and dust composition), se e IWolf et al] (1 1999b, Wolf & Hennind 
(120001) . IWolj (ho01[|2003h and I Schartmann et all (l2005fT 



1 MC3D (Monte Carlo 3D) has been te sted extensively again st other 
radiative transfer codes for 2D structures dPascucci et al .1120041) and we 
also performed a direct comparison for the special case of AGN dust tori 
with the simulations of Granato & Danes eli ll994h . one of the standard 
torus m odels for comparison, ca lculated with his grid based code (see 
Fig. 4 in lSchartmann et al.ll2005h . 



2.3. Resolution study 

In Fig. [2^, we show SEDs for our standard clumpy model (solid 
line, 5 ■ 10 6 monochromatic photon packages) and for the same 
model, but with twice as many photon packages (10 7 ) used for 
the simulation of the temperature distribution (dotted graphs, 
identical with solid lines). Despite slight differences in the tem- 
perature distributions of single grains, we find an almost iden- 
tical behaviour in the displayed SEDs, with differences smaller 
than the thickness of the lines. Maps at 12/im display the same 
distribution with slight problems along the projected torus axis, 
which are not visible in the single surface brightness distribu- 
tions and without any noticeable effect on the interferometric 
visibility distributions calculated from these maps. In Fig. |2j), 
the solid curves displays the SEDs for our high spatial resolu- 
tion standard model and the dotted lines refer to a model with a 
factor of roughly three less grid cells. Only very small deviations 
are visible at short wavelengths. Fig. [2] clearly shows that the re- 
sults and conclusions we draw from our simulations are neither 
affected by photon noise nor by too low spatial resolution. 



3. Analysis of our standard model 

In most of the SEDs discussed in this paper, only pure dust ree- 
mission SEDs are shown and an azimuthal viewing angle of 45° 
is used, if not stated otherwise. 



3.1. Temperature distribution 

In Fig. [3] the temperature distribution of all cells in all 9 and 
<p directions for the smallest silicate grain component is plotted 
a) for our clumpy standard model and b) for the corresponding 
continuous model. The red curves show the radial temperature 




r [pc] 

Fig. 3. Comparison between radial temperature distributions (for 
the smallest silicate grains) in all directions of the clumpy stan- 
dard model (panel a) with the temperature distribution of the cor- 
responding continuous model (panel b). The red curve indicates 
the temperature averaged over all angles 6 and <p. 
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profile, averaged over all (6, <p) directions. It is evident that the 
spread of temperature values for a given distance from the pri- 
mary radiation source is much larger for the clumpy models than 
for the continuous ones. Higher temperatures are possible even 
in parts of the torus further out, as dust free or optically thin 
lines of sight exist far out, depending on the distribution of sin- 
gle clumps. Therefore, a direct illumination of clouds is possible 
even at large radii. Concerning the continuous model, the scat- 
ter decreases significantly from 2pc outwards. Further in, the 8 
dependent radiation characteristic of the primary source causes 
greater scatter due to higher temperatures further away from the 
midplane. 

3.2. Viewing angle dependence 

Fig. [4^ shows the dependence of the spectral energy distribu- 
tions on the inclination of the torus. One can only see a clear 
distinction between lines of sight within the dust-free funnel (0° 
and 30° inclinations) and those within the wedge-shaped disk 
(60° a nd 90°). This was already reported bv lGranato & Danesd 
( 1 19941) . based on their continuous wedge models. In our case, 
it is caused by the relatively large volume filling fraction and 
the large clouds in the outer part of the torus. Therefore, only 
a weak dependence of the dust density distribution on the po- 
lar angle exists, which we chose for simplicity. In our previ- 
ous (2D) modelling (see ISchartmann et al.ll2.Q05b . we obtained 
the expected smooth transition in the polar direction. In Fig. |4j^, 
the azimuthal angle is varied for a constant inclination angle of 
60°. Nearly identical SEDs result, which is understandable when 
considering our large volume filling factors. The largest devia- 
tions appear at the shortest wavelengths, where the emission re- 
sults from the hottest parts of the torus, which are also the most 
centrally concentrated parts. Therefore, this wavelength range is 
most sensitive to changes of the optical depth along the direct 
line of sight towards the centre. 
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Fig. 4. Dependence of the SEDs on the viewing angle: a) dif- 
ferent inclination angles for a common azimuthal angle <p — 45°. 
Inclination angles shown are: 0° (solid line), 30° (dotted line), 
60° (dashed line), 90° (dashed-dotted line), b) Different az- 
imuthal angles <f> for a common inclination angle of i = 60°. 
Azimuthal angles shown are: 0° (solid line), 45° (dotted line) 
and 180° (dashed line). 
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Fig. 5. Inclination angle study of images of clumpy (first two 
rows) and continuous models (third and fourth row) with two 
different dust masses at 12 urn. Shown are the extreme cases with 
half of the mass of the standard model (first and third row) and 
eight times the mass (second and fourth row). For details of the 
mass study, see Sect. 14.31 The inclination angles i = 0°, 60°, 
90° are shown in different columns. The scaling is logarithmic 
with a range between the maximum value of all images and the 
10~ 6 th fraction of it (excluding the central point source). Labels 
are sizes in pc. 



The dependence on the inclination angle of images is shown 
in Fig. [5] (upper two rows). It is especially interesting that the 
different inclination angles look very similar, which was not the 
case for the continuous model (see lower two rows of Fig.0. 
There, the images at larger inclination angles are dominated by 
the boundaries of the disk, which are not so well defined in the 
clumpy case. In the zoomed-in images (Fig. [7] upper row), the 
basic features of our model are directly visible, as one can see 
the different illumination patterns of clouds: Clouds in the inner- 
most part are fully illuminated and, therefore, show bright inner 
rims and cold outer parts. Other clouds are partly hidden behind 
clouds further in and appear as bright spots only. 



3.3. Wavelength dependency 

Fig.|6]shows the wavelength dependency of our standard model. 
At short wavelengths, the hottest inner parts dominate the bright- 
ness distribution. Further out, a few more directly illuminated 
clumps are visible as bright spots. At the longest wavelengths, 
emission arises from clumps all over the torus, as colder dust 
emits strongest at these wavelengths. This dust is spread over 
a larger volume, due to the steeply decreasing temperature dis- 
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Fig. 6. Wavelength dependence of the surface brightness distri- 
butions: A = 4.6 um (upper row), A = 9.7 urn (second row), 
A = 12.0 um (third row) and A - 30.2 um (lower row). Within 
the rows, the inclination angle changes from face-on view (left- 
most panel) over 60° to 90° (rightmost panel). The images are 
given in logarithmic scaling with a range of values between the 
global maximum of all images and the 10~ 5 th fraction of it (cen- 
tral source excluded). Labels are in pc. 



tribution at small radii. Furthermore, the extinction curve has 
dropped by a large factor at these wavelengths and, therefore, 
the torus becomes optically thin and the whole range of cloud 
sizes is visible. 



4. Parameter variations 

4.1. Different realisations of the clumpy distribution 

As already discussed in lSchartmann et aTl (12005). SEDs of dust 
reemission depend strongly on the distribution of dust in the in- 
nermost region. Changing the random arrangements of clumps 
- as done in this section - therefore is expected to cause sig- 
nificant changes of the SEDs, especially for the case of a small 
number of clouds. The second important parameter is the opti- 
cal depth of the single clumps. The larger it is, the stronger is 
the dependence of the SEDs on the dust distribution in the in- 
nermost region. For the case of our modelling, the small number 
of clouds is expected to cause large differences in the observed 
SEDs. But this effect is partially compensated - in most of the 
simulations - by optically thin individual clumps, resulting in a 
more similar behaviour of the SEDs. 

Looking at the simulated SEDs and matching them with the 
images, the following results can be seen (compare to Fig. [7]): 



At 0° inclination angle (upper row), we observe nearly iden- 
tical SEDs. The dashed line (corresponding to the fourth col- 
umn) shows a slightly enhanced flux at short wavelengths, as a 
larger number of clouds are close to the central source. In the 
third column, the cloud number density in the central part is the 
lowest of the three examples. 

For the case of the middle row (i = 60°), the largest devi- 
ations are visible for the case of the dotted line (third column). 
Here, the silicate feature even appears in emission. This is visi- 
ble in the surface brightness distributions, as more directly illu- 
minated clouds are visible on unobscured lines of sight, result- 
ing in a brighter central region compared to the other two maps. 
At an inclination angle of 90° (third row), absorption along the 
line of sight increases drastically from the second to the fourth 
column, visible in a deepening of the silicate absorption feature 
and the darkening of the central region of the surface brightness 
distributions. 



4.2. Different volume filling factors 

Starting from the standard model with a volume filling factor of 
30% and 400 clumps within the whole model space, we halved 
it once by distributing only 160 clumps within the calculation 
domain and doubled it, for which 1500 clumps were needed due 
to the applied procedure of randomly distributing clumps. 

The resulting surface brightness distributions at A = 
12.0 um are shown in Fig.|8]with the three models given in differ- 
ent rows and for three different inclination angles: i = 0°, 60°, 
90°. In the case of the lowest volume filling factor, individual 
clouds are visible. The distribution of the surface brightness of 
these individual clouds reflects the temperature structure within 
single clumps. The directly illuminated clumps are hotter and, 
therefore, appear brighter. When adding more and more clumps 
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Fig. 8. Different volume filling factors: 15% (upper row), 30% 
(middle row), 60% (lower row) for A = 12.0 urn. From left to 
right, the inclination angle changes: i = 0°, 60°, 90°. The scaling 
is logarithmic with a range between the maximum value of all 
images and the 10~ 6 th fraction of it (excluding the central point 
source). Length scales are given in pc. 
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column of images, the dotted line to the second and the dashed line to the third column. Images are given in logarithmic color scale 
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(increasing the volume filling factor), the chance of directly illu- 
minating clumps further out decreases and at higher filling fac- 
tors it is only possible for clumps close to the funnel. This is 
clearly visible at higher inclination angles: the higher the vol- 
ume filling factor, the clearer the x-shaped feature appears, as 
only clumps within or close to the funnel can be directly illu- 
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Fig. 9. Different volume filling factors: 15% (dotted line), 30% 
(solid line), 60% (dashed line). Rows show SEDs for the three 
inclination angles: i = 0°, 60°, 90°. The viewing angle in in- 
direction is always 45°. 



minated. At a volume filling factor of 60%, the surface bright- 
ness distribution looks very similar to that of the corresponding 
continuous model (compare to Fig. |5). For large volume filling 
factors and close to edge-on, substructure is only visible from 
clouds in a viewing direction towards the dust-free cones. 

The corresponding SEDs are shown in Fig. [9] With increas- 
ing filling factor, more and more flu x at short wavelengths ap - 
pears for the face-on case, as seen in ISchartmann et aiTd2005l) . 
The shape of the clumpy model SEDs resemble the correspond- 
ing continuous model most (compare to Fig. [TUJ right column) 
for the highest volume filling factor. Concerning the silicate fea- 
ture, it increases slightly in emission as the amount of dust at the 
appropriate temperature increases as well (visible at the tran- 
sition from the lowest to the medium volume filling factor). 
The silicate absorption feature at higher inclinations strongly de- 
pends on the viewing angle (compare Fig. |9j) and c) especially 
for the model with the least number of clumps (dotted line). 
Thus, this study shows the validity of the simplification of using 
a smooth dust distribution in the case of very high torus volume 
filling factors, as was assumed in previous simulations. 

4.3. Dust mass study 

To study the dependence of the SEDs on the optical depth of 
the torus, we carried out a study with 0.5, 1, 2, 4 and 8 times 
the dust mass in the standard model. This leads to an opti- 
cal depth at 9.7 um within the equatorial plane, averaged over 
all angles of <p of \f^^f = 1, 2, 4, 8, 16. Single clumps 

then change from optically thin to optically thick (t^ 1 ™^ = 
0.19,0.38,0.76,1.52,3.04). 

The resulting behaviour of the SEDs is shown in Fig. [10] 
where it is also compared to the corresponding continuous mod- 
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Fig. 11. Close-up of the spectrum between 7 and 15 (xm of the 
face-on case of the dust mass study for the clumpy model in 
linear display. The lines are defined as in Fig.fTOl 



els. Concerning the silicate feature (in emission) for the face-on 
case (top row), we see a very similar behaviour of the SEDs 
of the clumpy and continuous model. Increasing the mass and 
with it the total optical depth leads to a flattening of the SED 
around the silicate feature, even more pronounced in the con- 
tinuous case. In addition to that, a slight shift of the maximum 
of the silicate feature towards longer wavelengths is visible for 
the case of the highest dust mass, apparent in the zoom-in of 
Fig.fTOl around the silicate feature (see Fig-ITTTi. This is due to the 
increasing underlying continuum towards longer wavelengths. 
Although the principal behaviour of the silicate feature is iden- 
tical for our clumpy and our continuous model, the reasons dif- 
fer: In the case of a continuous wedge-like torus, the inner, di- 
rectly illuminated walls are only visible through a small amount 
of dust. From step to step, the walls become more opaque and 
shield the directly illuminated inner rim better, decreasing the 
height of the silicate fe ature. This was not the ca se for our con- 
tinuous TTM- models in lSchartmann et aH (|2005j). With them, it 
was not possible to significantly reduce the silicate feature height 
within reasonable optical depth ranges. This was caused by the 
fully visible, directly illuminated inner funnel. Therefore, in the 
wedge-shaped continuous models, the reduction of the feature is 
an artefact caused by the unphysical, purely geometrically mo- 



tivated shape. Furthermore, it also involves very deep and so far 
unobserved silicate absorption features for the edge-on case (see 
lower right panel in Fig. [TOt. Concerning the clumpy model, the 
explanation for the flattening of the silicate feature in the edge- 
on case with increasing dust mass of the torus will be given in 
Sect. I5TT1 

The Wien branches show different behaviour when looking 
face-on. For the case of the clumpy torus model, increasing the 
optical depth means that the Wien branch moves to larger wave- 
lengths, as expected for the edge-on case. This is understandable 
when most of the directly illuminated surfaces of the clouds are 
then hidden behind other clouds, an argument which is not valid 
if the clouds are too optically thin in the inner part. 
For the edge-on case, we qualitatively obtain a comparable be- 
haviour as in the continuous case, because of the large number of 
clumps and the same optical depth within the equatorial plane. 
But a very important difference can be seen in the appearance 
of the silicate feature in absorption: when we want to have only 
very weak silicate emission features in the face-on case, a large 
optical depth is needed, resulting in an unphysically deep silicate 
feature in absorption in the edge-on case of the continuous mod- 
els, whereas the silicate feature remains moderate for many lines 
of sight for the clumpy model, where we see a large scatter for 
different random arrangements of clumps (compare to Fig. [7j- 

Concerning surface brightness distributions (see Fig.|5J, one 
can see that the objects appear smaller at mid-infrared wave- 
lengths for the case of higher dust masses: the larger the op- 
tical depth, the brighter the inner region and the dimmer the 
outer part. This is caused by a steepening of the radial dust tem- 
perature distribution with increasing mass of the objects, as the 
probability of photon absorption increases in the central region. 
Especially for the continuous case, the asymmetry at intermedi- 
ate inclination angles becomes visible for larger optical depths 
caused by extinction on the line of sight due to cold dust in the 
outer parts of the torus. 

4.4. Concentration of clumps in radial direction 

As already described in the model section ©, clump posi- 
tions are also chosen in accordance with the density distribu- 
tion of the corresponding continuous model. Therefore, chang- 



M. Schartmann et al.: Three-dimensional radiative transfer models of clumpy tori in Seyfert galaxies 



9 



-40 -20 



40 -40 -20 20 40 -40 -20 



0° ■ 


60° - 


90° - 

a 

Pi m w. ' 

a= 0.0 - 




; 


* 

a = -0.5 - 


• 




m 

a=-1.0 - 



-40 -20 



-40 -20 



Fig. 12. Surface brightness distributions for various slopes of 
the density distribution in the corresponding continuous model 
(Pcont r") leading to different concentrations of clumps in the 
radial direction in the clumpy model. The slopes are: a = 0.0 
(upper row), a = -0.5 (middle row), a = -1.0 (lower row). 
From left to right, the inclination angle increases from face-on 
to edge-on: i = 0°, 60°, 90°. Shown are images at 12 |im with a 
logarithmic color scale ranging from the maximum of all images 
to the 10~ 6 th part of it. Labels denote the distance to the centre 
in pc. 



The same behaviour is visible in the corresponding SEDs 
shown in Fig. [13] Decreasing the amount of dust in the centre 
near the heating source leads to decreasing flux at near-infrared 
wavelengths, whereas the flux at far-infrared wavelengths in- 
creases (reflecting the enhancement of dust in the outer part). 
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Fig. 13. SEDs for the clump concentration study. The varied 
slopes of the underlying density distribution are: a = 0.0 (dotted 
line), a = -0.5 (solid line), a = —1.0 (dashed line): a) face-on, 
b) edge-on. 



ing the slope of this radial density distribution, defined to be 
PcontC''; @> <P) = Po \t~z) > leads to a different concentration of 
clumps along radial rays. In this section, we vary the slope of the 
distribution a from a homogeneous dust distribution (a = 0.0) 
over a = -0.5 (our standard model) to a = -1.0. Decreasing a 
leads to an enhancement of the clump number density towards 
the central region. In order to keep the volume filling fraction 
at a constant level of 30%, we need to increase the number of 
clumps, as their size decreases towards the central region. All 
clumps possess the same optical depth. In order to have a con- 
stant mean optical depth in the midplane, the total dust mass has 
to be decreased. For an overview of the modified parameters see 
Table [2] The change of clump concentration can be seen directly 
from the simulated images at 12.0 [xm in Fig. [12] especially in the 
face-on case (first column). In the upper panel, single reemitting 
clumps are visible in the central region. This changes more and 
more to a continuous emission for the case of the highest cloud 
concentration in the centre due to multiple clumps along the line 
of sight and intersecting clumps. At higher inclination angles, 
the higher concentration leads to a sharper peak of the surface 
brightness. 



Table 2. Varied parameters of the clump concentration study. 



Parameter a = a = -0.5 a = -1 

No. of clumps 250 400 900~ 

Dust mass [M ] 22950 12562 6418 



4.5. Dependence on the clump size distribution 

In our clumpy standard model, a radially changing clump size 
proportional to the radial distance to the centre was chosen. In 
this section, we test the effects of decreasing the slope /3 of the 

radial size distribution a c i ump = ciq (yjjjr J 8 of the clumps. This is 
done in a way that the volume filling fraction as well as the op- 
tical depth in the midplane, averaged over all azimuthal angles 
<p, remain constant. It is achieved by changing the proportional- 
isation constant of the clump size distribution ao and the total 
dust mass of the torus. Doing this results in very well resolved 
clumps in the inner part. Beyond a distance of approximately 
25 pc, the number of grid cells per clump drops below the value 
of our standard model. 

Surface brightness distributions for the extreme case of a 
constant clump size are shown in Fig. [14] Due to the large clump 
radius of 5pc even in the central region and that fractions of 
clumps at the model space boarder are prohibited, it leads to a 
density distribution with a quite large, unevenly shaped central 
cavity, as can be seen in the face-on view (left panel of Fig. PHI , 
The inner rim is given by only a few intersecting clumps, instead 
of the otherwise defined spherical central cavity. Therefore, in 
the edge-on case, the surface brightness distribution shows an in- 
ner boundary, which is bent towards the centre (convex shaped). 
In these models, due to the large clump size in the inner region, 
many clumps intersect, producing a nearly continuous dust dis- 
tribution at the inner boundaries, which lets the - typical for con- 
tinuous models - x-shaped structure appear again. For the same 
reason, the extinction band due to the | cos(#)|-radiation char- 
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acteristic is visible in the edge-on view. Especially at the 90° 
inclination angle, single clumps are directly visible (above and 
below the centre). In these cases, their shading directly shows the 
illumination pattern due to the primary source (accretion disk), 
emission from other clumps and extinction from the foreground 
dust distributions. 

The corresponding SEDs (Fig. fT3T > mainly reflect the in- 
crease of the inner cavity and, therefore, the lack of flux at short 
wavelengths. The convex shape of that region causes a larger di- 
rectly illuminated area at the funnel walls and, therefore, slightly 
strengthens the silicate emission feature in a face-on view. A dif- 
ferent appearance (emission/absorption) of the 10 urn feature at 
i = 60° (middle panel) is seen. This is due to the lower number 
density of clumps in the inner part, enforced by the restriction of 
having only whole clumps within the model space. 

A dust mass study for the case of the large, constant diameter 
clump model (J3 = 0) reveals the same behaviour as discussed 
in Sect. 14. 31 when looking edge-on onto the torus. However, the 
face-on case differs: only the relative height of the silicate feature 
c hanges slightly. This was already observed in our 7TM-models 
in lSchartmann et all d2005l) and is due to the now inwardly bent 
inner walls of the funnel (see Fig.[14j right panel), caused here 
by the very large and spherical clumps in the innermost torus 
region. 

5. Discussion 

5.1. Explanation for the reduction of the silicate feature 

The results shown in the subsections before can be explained 
with the following m odel, which was partially discussed by 
iNenkova et alj d2002l) . It is illustrated in Fig. [16] where yellow 
denotes directly illuminated clump surfaces. Many of the ex- 
plained features can also be seen in the zoomed-in versions of 
the surface brightness distributions (Fig.|7]i for the face-on case 
(upper row). 

As already pointed out in lSchartmann et ai1d2005l) . the SEDs 
of dust tori in the mid-infrared wavelength range are mainly de- 
termined by the inner few parsecs of the toroidal dust distribu- 
tion. In each of the central clouds of the clumpy model, the dust 
temperature drops from the inner directly illuminated edge to- 
wards the cloud's outer surface. 

With an inclination angle close to i = 90°, we expect - for 
realistic volume filling factors - comparable behaviour of the 
SED as in the continuous case. The silicate feature has a smaller 
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Fig. 14. Images at 12 um and inclination angles i = 0° and 90° 
for a clumpy distribution with a constant clump size, indepen- 
dent of the radial position. The scaling is logarithmic with a 
range between the maximum value of all images and the 10~ 6 th 
fraction of it (excluding the central point source). 
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Fig. 15. Dependence on the clump size for different inclination 
angles (rows: 0°, 60°, 90°). The solid line corresponds to our 
standard model with /3 = \. For the dashed line, clumps have 
equal size (a c iump = 5 pc, f3 = 0), independent of the radial po- 
sition and the dotted line corresponds to an intermediate model 
with/3 = 0.5. 

depth, as discussed in Sect. 14. 31 But the situation changes with 
decreasing inclination angle. Here, one has to distinguish be- 
tween different cases: 

1. With a relatively high volume filling factor and a not too 
small extension of the clouds in the central region of the 
torus, it is likely that the directly illuminated part of most 
of the clouds is hidden from direct view by other clouds. 




Fig. 16. Sketch of our clumpy torus model. Indicated in yellow 
are directly illuminated surfaces of the clumps, i is the inclina- 
tion angle, open is the half opening angle of the torus. 
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Therefore, the directly illuminated surface area is reduced 
compared to the corresponding continuous model. As this 
area is responsible for the emission fraction of the silicate 
feature within the SED, it shows less silicate emission. In 
order to produce such a shadowing effect, clouds have to 
possess a large enough optical depth, which means that they 
have to be either small or massive (r clump oc ra c i ump a^ ump , 
where m c i ump is the mass of the clump). 

2. For the case of a too small optical depth of the clumps in the 
innermost region, we expect the silicate feature to appear in 
strong emission. 

3. Another possibility of producing silicate features in emission 
is when the model possesses a small number of clumps in the 
inner part, making the shadowing effect inefficient. 

4. A further effect on the silicate feature strength arises from the 
grain dependent sublimation implemented in our models. As 
graphite grains possess a higher sublimation temperature as 
silicate grains, they are able to partially shelter the silicate 
grains from direct irradiation. 

Thus, the strength of the silicate feature is mainly deter- 
mined by the distribution, size and optical depth of the clouds 
in the direct vicinity of the sublimation surface of the dust. We 
will see in the next sec t ion th at thi s finding well explains the 
fact th at iNenkova et al.l d2002l) and iDullemond & van Bemmell 
d2005l) come to different conclusions concerning the reduction 
of the strength of the silicate feature due to dumpiness, as the 
distribution (and size) of their clouds differ. 

5.2. Comparison with other torus models 

The results of lNenkova et al.1 d2002l) - the pioneering work in the 
field of clumpy tori - are broadly consistent with the explanations 
giv en in Sect.l5.1lof this paper. 

IDullemond & van Bemmell d2005l) model 2D clumps in the 
form of rings with a two-dimensional radiative transfer code. In 
contradiction to all other simulations, no systematic reduction 
of the silicate feature due to dumpiness is found. The reason for 
this is understandable with the explanations given in Sect. 15.11 as 
their model features a small clump number density in the central 
region and shadowing effects are rather small. Therefore, they 
find both strengthening of the silicate feature and reduction, de- 
pending on the random ring distribution. 

A comparison of our clumpy standard model a nd two other 
random cloud distributions with simulations by IHonig et"aT] 
d2006l) is shown in Fig. [17] They follow a different, multi- 
step approach: 2D radiative transfer calculations of individual 
clouds at different positions and with various illumination pat- 
terns within the torus are carried out. In a second step, the SED of 
the total system is calculated. The cloud distribution and parame- 
ters such as optical depth or size arise from an accr etion scenario 
of self-gravitating clouds close to the shear limit dVollmer et al.l 
l2004t iBeckert & DuschHl2004l) . The advantage of this approach 
is that resolution problems can be overcome easily, as only 2D 
real radiative transfer calculations of single clumps are needed. 
Characteristics for their modelling are small cloud sizes with 
very high optical depths in the inner part of the torus and a large 
number of clumps. For comparison, a cloud at the sublimation 
radius of their model has a radial size of R c i ou d = 0.02 pc with 
an optical depth of Tq 1 ™^ 



sity there reduces the silicate feature significantly by shadowing 
with respect to their single clump calculations. Their finding that 
the silicate emission feature can be reduced further by increas- 
ing the number density of clumps in the innermost part perfectly 
fits our explanation presented in Sect. 15.11 Deviations between 
the two approaches (see Fig.[T7ll are mainly due to the approxi- 
mately eight times larger primary luminosit y and the larger op- 
tical depth, at least in the midplane of the IHonig et~aT] d2006l) 
modelling compared to our standard model. Furthermore, in our 
simulations only dust reemission SEDs are shown. This leads to 
relatively higher fluxes at short wavelengths compared to long 
wavelengths for the i = 30° case (Fig. [TTh ) and to more ex- 
tinction within the midplane and, therefore, a shift of the Wien 
branch towards longer wavelengths in the edge-on case (lower 
panel). 

6. MIDI interferometry 

Even with the largest single-dish mid-infrared telescopes, it 
is impossible to directly resolve the dust torus of the nearest 
Seyfert galaxies. Therefore, interferometric measurements are 
needed. Recently, iJaffe et ail (120041) succeeded for the first time 
to resolve the dusty structure around an AGN in the mid-infrared 
wavelength range. In this case, they probed the active nucleus of 
the nearby Seyfert 2 galaxy NGC 1068 with the help of the MID - 
infrared interferometric Instrument (MIDI. ILeinert et al1l2003l) . 
It is located at the European Southern Observatory's (ESO's) 
Very Large Telescope Interferometer ( VLTI) laboratory on Cerro 
Paranal in Chile. Its main objective is the coherent combination 
of the beams of two 8.2 m diameter Unit Telescopes (UTs) in 
order to obtain structural properties of the observed objects at 
high angular resolution. A spatial resolution of up to A/ (2 B) » 
10 mas at a wavelength of A = 10 urn can be obtained for the 
largest possible separation of two Unit Telescopes of B « 120 m. 
Operating in the N-band (8 - 13.5 urn), it is perfectly suited to 
detect thermal emission of dust in the innermost parts of nearby 



10 J 



JS r 



1 cr 



icr 







i = 30° - 
\ 




^ 


1 = 90° - 



10 



100 



250. In our standard model, clouds 
at the sublimation radius are four times larger and possess an 
optical depth of only Tq 1 ™^ ~ 3. The large optical depth in 
the innermost part in combination with the large number den- 



Fig. 17. Comparison of our clumpy standard model and two 
other random realisations of the clump distribution (blue lines) 
with simulat ions done by S. F . Honig (private communication, 
described in IHonig et al. 2006), shown by the yellow lines, for 
10 different random realisations of their model. The latter are 
scaled with a factor of 2.2 in order to give rough agreement be- 
tween the two models (see text for further explanation). 
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Fig. 18. Visibilities of our continuous standard model at a wave- 
length of 12 (xm plotted against the projected baseline length. 
Colour of the visibility distributions refers to different position 
angles of the projected baseline w. r. t. the torus axis. Each panel 
shows a different inclination angle, as indicated in the upper right 
corner. 
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Fig. 19. Visibilities of our clumpy standard model at a wave- 
length of 12 u.m plotted against the projected baseline length. 
Colour of the visibility distributions refers to different position 
angles of the baseline w. r. t. the torus axis. Each panel shows a 
different inclination angle, as indicated in the upper right corner. 



Seyfert galaxies. MIDI is designed as a classical Michelson in- 
terferometer. Being a two-element beam combining instrument, 
it measures so-called visibility amplitudes. Visibility is defined 
as the ratio between the correlated flux and the total flux. Its in- 
terpretation is not straightforward, since no direct image can be 
reconstructed. Therefore, a model has to be assumed, which can 
then be compared to the visibility data. MIDI works in dispersed 
mode, which means that visibilities for the whole wavelength 
range are derived. The dust emission is probed depending on the 
orientation of the projected baseline. Point-like objects result in 
a visibility of one, as the correlated flux equals the total flux. 
The more extended the object, the lower the visibility. With the 
help of a density distribution, surface brightness distributions in 
the mid-infrared can be calculated by applying a radiative trans- 
fer code. A Fourier transform of the brightness distribution then 
yields the visibility information, depending on the baseline ori- 
entation and length within the so-called U-V-plane (or Fourier- 
plane). 

The main goal of the following analysis is to investigate 
whether MIDI can distinguish between clumpy and continuous 
torus models of the kind presented above. Furthermore, we try to 
derive characteristic features of the respective models and show 
a comparison to data obtained for the Circinus galaxy. 

6. 1 . Model visibilities 

In Fig. [18] calculated visibilities for four inclinations of our con- 
tinuous standard model at a wavelength of A — 12 urn are shown. 
Various orientations of the projected baseline are colour coded 
(the given position angle is counted anti-clockwise from the pro- 
jected torus axis). Due to the axisymmetric setup, all lines coin- 
cide for the face-on case. For all other inclination angles, visibili- 



ties decrease until a position angle of 90° is reached and increase 
symmetrically again. This means that the torus appears elon- 
gated perpendicular to the torus axis at this wavelength. Fig. [T9l 
shows the same study, but for the corresponding clumpy model. 
The basic behaviour is the same, but the visibilities show fine 
structure and the scatter is much greater, especially visible in 
the comparison of the i = 0° cases. Furthermore, while all of 
the curves of the continuous model monotonically decrease with 
baseline length, we see rising and falling values with increasing 
baseline length for the same position angle in the clumpy case. 
In addition, for the continuous models, curves do not intersect, 
in contrast to our clumpy models. However, to detect such fine 
structure in observed MIDI data, a very high accuracy in the vis- 
ibility measurements of the order of <x v ~ 0.02 and a very dense 
sampling is required. 

In Fig. [20] the wavelength dependence of the visibilities is 
shown. The first two panels represent the case of the clumpy 
standard model and the third and fourth the continuous model. 
Each of the two panels of the respective model visualises a dif- 
ferent position angle (counted anti-clockwise from the projected 
torus axis). An inclination angle of 90° is used in all panels. 
Three different wavelengths are colour coded: 8.2 urn at the be- 
ginning of the MIDI-range (black dotted line), 9.8 urn within the 
silicate feature (blue) and 12.6 urn at the end of the MIDI wave- 
length range (yellow), outside the silicate feature. While the con- 
tinuous model results in smooth curves (see also Fig.[T8l, much 
fine structure is visible for the case of the clumpy model. The 
differences between the displayed wavelengths relative to the 
longest wavelength are smaller for the clumpy models than in 
the continuous case. 

Fig. [2T] shows visibilities for our clumpy standard model 
at 12 urn, plotted against the position angle (counter-clockwise 
from the projected torus axis). Baselines are colour coded be- 
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Fig. 20. Visibilities of our clumpy (first two panels) and continuous (last two panels) standard model at different wavelengths 
(colour coded), plotted against the projected baseline length for the two position angles 0° (in torus axis direction) and 90° (along 
the midplane) for an edge-on view onto the torus (i = 90°). 
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Fig. 21. Visibilities of our clumpy standard model at different 
inclination angles (as annotated in the upper right corner) plotted 
against the position angle for various projected baseline lengths 
(colour coded) and a wavelength of 12 (im. 



tween 20 m and 100 m in steps of 20 m. A longer baseline means 
that structures are better resolved, leading to decreasing visibili- 
ties. For the case of inclination angles close to edge-on, the vis- 
ibility distribution changes from more or less flat to a character- 
istic oscillating distribution at longer baselines (from 60 m on- 
wards) with minima around 100° and 300°. This means that our 
torus model seems to be more elongated within the equatorial 
plane and has the smallest width along the projected torus axis 
at this wavelength. But this only applies for the innermost part; 
the torus as a whole looks approximately spherically symmetric. 
At small inclination angles no such favoured size distribution is 
visible. 

6.2. Comparison with MIDI-data for the Circinus galaxy 

Unfortunately, a fitting procedure involving a large parameter 
study is not possible with our current model, due to the very 



long computational times of the order of 30 to 40 hours per in- 
clination angle (including calculation of the temperature distri- 
bution, the SED and surface brightness distribution). Therefore, 
we applied the following procedure: From our experience with 
modelling the SED of the Circinus galaxy w ith our previously 
used continuous Turbulent Torus Models (see lSchartmann et al.l 
120051) . we adopt the size of the object used there. Furthermore, 
we tried to stay as close to our clumpy standard model as possi- 
ble (for the parameters of the clumpy standard model, compare 
to Table[TJ and copied the parameters a, j3 and oq. The rest of 
the parameters were changed, in order to obtain the best pos- 
sible adaptation to the data, within the investigated parameter 
range. The comparison of our current clumpy Circinus model 
as described above and in Table [3] (yellow stars) to interferomet- 
ric observations with MIDI dTristram et al.ll20 07) of the Circinus 
galaxy (black) is shown in Fig. [22] In contrast to the presenta- 
tion of continuous visibility curves above, single measurements 
of combinations of various baseline lengths and position angles 
are displayed in this plot. Position angle now refers to the an- 
gle on the sky measured from north in a counter-clockwise di- 
rection. The rotation axis of our simulated torus has a position 
angle of approximately -45° according to this definition. The 
black numbers denote the length of the projected baseline (given 
in m) of the corresponding data point. From the approximate 
correspondence of the model values with the data, one can see 
that the size of the emitting region at the two wavelengths is re- 
produced quite well. Most of the local extrema of the curve can 
be reproduced for the case of 9.1 \im. Larger deviations are vis- 
ible for A = 12.0 ^m. The good adaptation partly is due to the 
changes in baseline length. Longer baselines naturally result in 
smaller visibilities, as we are probing smaller and smaller struc- 
tures (see also Fig.[l9]l. Greater visibilities for shorter or equal 



Table 3. Circinus model parameters: For an explanation of the 
parameters see Sect.[2]and Table [T] Mbh is the mass of the cen- 
tral black hole (from lGreenhill etaT]|2003l) and Ldisk/^edd is the 
Eddington luminosity ratio resulting for the assumed luminosity 
of the central source. 
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baselines and similar position angle, therefore, have to be due 
to those curves in Fig. [19] with increasing visibility with base- 
line or a very inhomogeneous distribution of dust with position 
angle. Both can be interpreted as signs of dumpiness. The SED 
of the same Circinus model is plotted over current high resolu- 
tion data in Fig. [23] The NIR (near-infrared) data points were ob- 
tained with the NACO camera at th e VLT and correcte d for fore- 
ground extinction by Ay — 6 mag dPrieto et alj |2004 : ). Different 
symbols refer to various aperture sizes (se e figure caption). The 
thick green line shows the MIDI spectrum (iTristram et al.|[2007h 
and the black line is our Circinus model as discussed above for 
an aperture of 0.4" in radius; the yellow line denotes the same 
model, but calculated for the whole simulated model space. Both 
modelled SEDs include the direct radiation of the central source 
(calculated with real Monte Carlo radiative transfer), which in 
these examples dominates over dust reemission for the small 
wavelength part from about 2 to 3 urn downwards and shows 
some noise, due to the low photon packet num bers used. In con- 
tradict ion to our continuous Circinus model in lSchartmann et all 
(120051) . enough nuclear radiation can be observed in order to ex- 
plain the turnover of the SED at small wavelengths and we do 
not need to assume scattering by material (dust and electrons) 
within the torus funnel. As can be seen from these figures, our 
model is able to qualitatively explain the SED as well as the visi- 
bility information. However, as we are not able to investigate the 
whole parameter range of our models, we cannot exclude that a 
different parameter set can describe the data equ ally well. This 
degen eracy problem was already pointed out by iGalliano et al.l 
d2003l) for the case of SED fitting. Adding new dumpiness pa- 
rameters will even strengthen this degeneracy. On the other hand, 
adding more data such as more visibility information will place 
more constraints and will weaken this problem. 

7. Conclusions 

In this paper, we implemented a new clumpy torus model in three 
dimensions. For computational reasons, a wedge-like shaped 
disk is used. In the discussion of our results, we place special 
emphasis on the comparison with continuous models and their 
differentiation using interferometric observations, such as with 
MIDI 

In lSchartmann et alj d2005l) . we had found that the SEDs of 
AGN tori in the mid-infrared wavelength range are mainly de- 
termined by the innermost part of the torus. With the presented 
clumpy torus models, this claim can be further strengthened. 
According to the new simulations, the silicate feature strength 
is mainly determined by the number density and distribution, as 
well as the optical depth and size of the clumps in the inner re- 
gion. With a sufficiently high optical depth of the clouds in the 
inner part, shadowing effects become important, which hide the 
illuminated cloud surfaces from direct view and, thereby, reduce 
the silicate feature in emission. At the same time, enough lines 
of sight with low optical depth remain so that only weak absorp- 
tion features result for the edge-on case. Continuous models with 
special and unrealistic morphologies (like the wedge-shaped tori 
used here) are also able to weaken the silicate emission fea- 
ture for the face-on view when applying an anisotropic radiation 
characteristic, but fail to simultaneously account for moderate 
absorption features, when looking edge-on to the torus. 

Due to the large clumps in our model, appreciable scatter 
in SEDs for different random realisations of the torus are ex- 
pected. A contrary effect is caused by the small optical depth 
of the single clumps and also of many dust-free lines of sight 
towards the centre. Direct comparison between calculated inter- 



ferometric visibilities for clumpy and the corresponding contin- 
uous models show that clumpy models naturally possess more 
fine structure, which can partly be resolved by MIDI. 

We also showed that these kinds of models are able to quali- 
tatively describe the available interferometric visibility and high 
resolution spectroscopic data of the Circinus galaxy at the same 
time. Currently, it is one of the best studied S eyfert galaxies 
in ter ms of mid-infrared visibility measurements dTristram et al.l 
|2007|) . The decreasing slope of the SED at short wavelengths can 
be described with our clumpy mod el, whereas it was at odds with 
the continuous model described in lSchartmann et ail d2005l) . 
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Fig. 22. Comparison of model visibilities (yellow stars and 
lines) for an azimuthal viewing angle of <p — 225° with 
MIDI observations for two different wavelengths. The base- 
line length for all data points is given above the data. Data 
courtesy of lTristram et alj d2007l) . 
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Fig. 23. Comparison of model SEDs with data for the Circinus 
galaxy. Different symbols refer to various aperture radii: blue 
stars - 0.38" (NACO), red rectangle - 0.1" (HST/NICMOS) 
an d black triangle - 1.0" (ESO/TIMMI2). Data compilation 
by iPrieto et al.l 12004 The thick gre en line shows the total 
MIDI spectrum dTristram et al.|[2Q07h . Our model (see model 
parameters in Table|3j is calculated for an aperture radius of 
0.4" (black line) and the total model space (yellow line). 



